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Abstract: We formulate the high temperature expansion in supersymmetric matrix quan- 
tum mechanics with 4, 8 and 16 supercharges. The models can be obtained by dimen- 
sionally reducing J\f = \ U(A^) super Yang- Mills theory in D = 4, 6, 10 to 1 dimension, 
respectively. While the non-zero frequency modes become weakly coupled at high temper- 
ature, the zero modes remain strongly coupled. We find, however, that the integration over 
the zero modes that remains after integrating out all the non-zero modes perturbatively, 
reduces to the evaluation of connected Green's functions in the bosonic IKKT model. We 
perform Monte Carlo simulation to compute these Green's functions, which are then used 
to obtain the coefficients of the high temperature expansion for various quantities up to 
the next-leading order. Our results nicely reproduce the asymptotic behaviors of the re- 
cent simulation results at finite temperature. In particular, the fermionic matrices, which 
decouple at the leading order, give rise to substantial effects at the next-leading order, re- 
flecting finite temperature behaviors qualitatively different from the corresponding models 
without fermions. 
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1. Introduction 

Recently large- A'^ gauge theories are playing increasingly important roles in string theory. 
One of the crucial discoveries was that U(A^) gauge theory appears as a low energy effective 
theory [1] of a stack of N D-branes [2] in string theory. This led to various interesting conjec- 
tures. For instance, large- iV gauge theories obtained by dimensionally reducing lOd U(A'') 
super Yang-Mills theory to 0, 1, 2 dimensions are conjectured to provide non-perturbative 
formulations of superstring/M theories [3-5]. 

Another type of conjectures asserts the duality between strongly coupled large- AT gauge 
theory and weakly coupled supergravity. In the AdS/CFT correspondence [6], for instance, 
it is conjectured that 4d = 4 U(A^) super Yang-Mills theory is dual to the type IIB 
supergravity on AdSs x S^. This duality is generalized to the finite temperature setup [7] and 
to non-conformal gauge theories [8]. Motivated by such dualities, large-iV gauge theories 
at finite temperature [9-12] have been studied intensively. 

Monte Carlo simulation is expected to be a powerful approach to explore the phase 
diagram of large- gauge theories. Indeed there was a remarkable progress in this direction 
recently. Supersymmetric matrix quantum mechanics have been studied by Monte Carlo 
simulation for the first time [13-15]. Ref. [15], in particular, deals with the model with 16 
supercharges, which may be viewed^ as the low energy effective theory of N DO-branes in 
the type IIA superstring theory [1]. The Monte Carlo results confirmed the gauge/gravity 
duality from first principles. Unlike in the bosonic model [10, 16, 17], no phase transi- 
tion was observed at finite temperature, which is consistent with the prediction based 

^The model is formally the same as the non-perturbative formulation [3] of the M theory although the 
large- iV limit should be taken in a different way. 
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on the gauge/gravity duality [11, 18]. Moreover, the internal energy at low temperature 
agreed with that of the non-extremal black hole described by the dual geometry. This 
implies in particular that the Bekenstein-Hawking entropy of the black hole is given a 
microscopic origin in terms of the open strings attached to the constituent N DO-branes. 
Unlike in Strominger-Vafa's result [19] for the extremal black hole, which relied on the 
supersymmetric non-renormalization theorem, the agreement has been found by studying 
the strongly coupled dynamics of the DO-brane effective theory directly. See refs. [8, 20] 
for earlier works, which connect the supersymmetric matrix quantum mechanics at finite 
temperature to the black-hole physics through the gauge/gravity duality. 

In this paper we formulate the high temperature expansion in the supersymmetric 
matrix quantum mechanics. While the low temperature behavior of the theory describes the 
classical black hole, the high temperature behavior is expected to describe hot strings [21]. 
We study the models with 4, 8, 16 supercharges that can be obtained by dimensionally 
reducing M = 1 U(A^) super Yang-Mills theory in L> = 4, 6, 10 to 1 dimension. The high 
temperature limit of the D = 10 case [3] for A/" = 2 has been studied in ref. [22] . As observed 
there and also in refs. [10,23], only the bosonic zero modes survive at the leading order, and 
their dynamics arc governed by the bosonic part of the IKKT [4] matrix model. In order to 
see the effects of the fermions, we proceed to the next-leading order. After integrating out 
the weakly-coupled non-zero frequency modes pcrturbativcly, we find that the remaining 
integration over the zero modes reduces to the evaluation of connected Green's functions in 
the bosonic IKKT model. This can be done by Monte Carlo simulation with much less effort 
than simulating the supersymmetric models at finite temperature directly. In particular, 
we are able to make a reliable large-A^ extrapolation using the data for N up to 32. We 
calculate the internal energy, the Polyakov line, and the extent of the eigenvalue distribution 
explicitly for D = 4, 6, 10. Our results nicely reproduce the asymptotic behaviors of the 
recent Monte Carlo data obtained for both supersymmetric models and bosonic models 
at finite temperature. The different properties of the two classes of models are clearly 
reflected in the next-leading order terms. 

The rest of this paper is organized as follows. In section 2 we define the model and the 
observables we study in this paper. In sections 3 and 4 we present the calculations at the 
leading order and at the next-leading order, respectively. Section 5 is devoted to a summary 
and discussions. In Appendix A we derive a formula, which is used to calculate the internal 
energy. In Appendix B we present the form of Green's functions used to evaluate them 
efficiently in actual Monte Carlo simulation. 

2. The models 

The models we study in this paper are defined by the action 

^^7 {^(^*^^)' + ^V'aA^a - \[Xi,Xj]'' - ^Mli)ap[Xi,^p]^ , (2.1) 

where Dt = dt — i[A{t), ■ ] represents the covariant derivative. The bosonic matrices 
A{t), Xi{t) (i = 1, 2, • • • ,d) and the fermionic matrices tpait) (« = 1) 2, • • ■ ,p) are N x N 
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Hermitian matrices, where p = 4, 8, 16 for d = 3, 5, 9, respectively. The models can be 
obtained formally by dimensionally reducing M = 1 super Yang-Mills theory in D = d + 1 
dimensions to one dimension, and they can be viewed as a Id gauge theory, where A{t), 
Xi{t) and tpait) are the gauge field, adjoint scalars and spinors, respectively. The p x p 
symmetric matrices ji obey the Euclidean Clifford algebra {7j,7j} = ^Sij. We impose 
periodic and anti-periodic boundary conditions on the bosonic and fermionic matrices, 
respectively. The extent f3 in the Euclidean time direction t then corresponds to the 
inverse temperature /? = T~^. 

The action is invariant under the shifts 



A{t) ^ A{t) + a{t)l , Xi{t) ^ Xi{t) + Xil , (2.2) 

where a{t) is an arbitrary periodic function and Xj is an arbitrary constant. In order to 
remove the corresponding decoupled modes, we impose the conditions 

/•/? 

tr^(i) = 0, / d.ttTXi{t)=0 {i = l,2,--- ,d) . (2.3) 

Jo 

The 't Hooft large- A'^ limit corresponds to sending A?^ to oo with the 't Hooft coupling 
constant A = g^N fixed. Since the coupling constant g can be absorbed by rescaling the 
matrices and t appropriately, we can set A to unity without loss of generality. This implies 
that we replace the prefactor ^ in the action (2.1) by TV in what follows. 

We define the extent of the eigenvalue distribution and the Polyakov line as 

= ^ f^dtti(Xi{t)y , (2.4) 

/3 



NpJo 

P = ^tiVexp (^i J dtA{t)^ , (2.5) 



where the symbol "7^ exp" represents the path-ordered exponential. 

As a fundamental quantity in thermodynamics, the free energy J- = —^lnZ((3) is 
defined in terms of the partition function given in the present model as 

Z{(3) = J [VA]fs[VXUV^P](,e-^^^^ , (2.6) 

where the suffix of the measure [ ■ J/j represents the period of the field to be path-integrated. 

However, the evaluation of the partition function Z{p) is not straightforward in Monte 
Carlo simulation, which wc use for the integration over the zero modes. We therefore 
study the internal energy defined by 

which has equivalent information as the free energy, given the boundary condition = 
E at T = 0. Note also that the internal energy at T = provides the ground state 
energy of the quantum mechanical system, which should vanish unless the supersymmetry 
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is spontaneously broken. In Appendix A we show that the internal energy E can be 
expressed as 

§2 = (^b) + (^f) , (2.8) 
where the operators ^5 a^nd are defined by 

3 1 



= {i^a{li)ap[Xi,^p\) . (2.10) 



The symbol ( ■ ) represents the expectation value with respect to the model (2.1). 
Let us take the static gauge dtA{t) = 0. Correspondingly we add the ghost term 



S^^ = N j^dtiT(dtc{t)Dtc{t)) (2.11) 



to the action, where c, c are N x N matrices representing the ghosts. We make a Fourier 
expansion of the fields as 

Xi{t) = Y,K exp(ina;i) , Va(0 = XI exp(ira;i) , (2.12) 

n r 

c{t) = Cn exp^incot) , c{t) = Cn exp{—inLot) , (2-13) 

where u = ^ represents the unit of Matsubara frequencies, and the indices n and r take 
integers and half-integers, respectively, due to the imposed boundary conditions. In terms 
of the Fourier modes, the gauge-fixed action is written as 

S = So + S'kin + ^int , (2-14) 

So ^ -A^/?tr|i([AX^])' + ^([X^,X^])'} , (2.15) 



5kin = Nf3 trj i X(na;)2xi„X; + J2('^^)^^nCn + ^ ira;^-.Vr| , (2.16) 
5int = -NP trj J2 nujXl^[A, X^] + ^ na;Cn[yl, c„] + ^ ^ Vr] (2.17) 

r,s n^O npq 

where the symbol implies that the m = n=p=0 term is excluded. 
3. Leading order calculation 

In this section wc consider the high temperature limit [10,22,23], which corresponds to the 
leading order calculation at high T. From (2.14), one can easily see that all the non-zero 
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modes decouple, and one is left with the zero modes governed by the action (2.15). By 
rescaling the zero modes as 

A = T-^'^Xi (^ = 1, 2, • • • , d) , Ad = T'^'^A , (3.1) 

where we recall that D = d + 1, the zero-mode action can be brought into the form 

So = ^Ntv {F^^f , F^, = -i [4, A,] . (3.2) 

Here and henceforth, the Greek indices /x, u are assumed to run over 1, 2, • • • ,D. The 
dimensionally reduced (DR) model (3.2) is nothing but the bosonic part of the IKKT type 
matrix model [24-26] . The leading behavior of the observables at high temperature can be 
obtained as 

(i?^>:.TV2^1tr(i,)^^^^ = XiTV2, (3.3) 
(P)-l4T-3/^(^tr(i.)^)^^^l-lx.T-3/^ (3.4) 
^•Ec.'-T(^tr{F,A =^X2T, (3.5) 



7V2 4 \7V ^ /^^ 4 

where ( • )dr represents the expectation value with respect to the DR model, and the 
coefficients xi and X2 are given as 



1 



1 ^..^f^ ^2\ dC2 / 1 , ^2\ fj -,^ 1 



We note that the expectation values appearing here are standard quantities calculated in 
the DR model (3.2) by various methods [25,26]. In particular, the quantity in eq. (3.7) 
can be obtained exactly by simply rescaling the dynamical variables or by writing down 
the Schwinger-Dyson equation [25]. 

4. Next-leading order calculation 

As is clear from the previous section, the leading order results arc insensitive to the existence 
of fermions. In order to sec their effects, we need to proceed to the next-leading order 
calculation, which involves the integration over the non-zero modes. 
For that purpose, let us rescale the non-zero modes as 

Xl = Xl, A = A, Cn = Cn, ^ = ^^'^ In , (4.1) 

where n 7^ 0, so that the kinetic terms take the canonical form 

5Mn = A^tr|^^(27rn)2xi„X;-|-^(27rn)2g„c„ + i^27rir^_,^,| . (4.2) 
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Then the propagators are given by 

{'4^ar)ab{i^l3s)cdjj = 2T^irN '^"/^'^^ -«'^«<^'^^c , (4.4) 

where the symbol (( • )) represents integrating only the non-zero modes using the quadratic 
terms (4.2). The interaction terms are given by S'int = — Yli=i where 



Vi = V^TV J] ntr {xI^[Ad, K\) , V2 = V^Nj2ntj: (dn[AD, c„]) , 

V3 = ^ViAT^tr (lP_r[AD,4>r]) , V4 = ^V^iV^tr (^_^7i[ii, V^r-]) , 
r r 

V6^\eN^t. {[A, XL J [A, Xi] + [A, J [XI, i,]) . (4.6) 



Wc have introduced the expansion parameter^ e = and omitted terms, which are 

irrelevant to the calculations at the next-leading order. 

First we calculate the extent of the eigenvalue distribution (2.4), which can be decom- 
posed as 

= Itr {Xlf + ^ E (^n^-n) ■ (4-7) 

Let us consider the first term. The leading order contribution is given by (3.3). At the 
next-leading order, we use the vertices (4.6) and integrate over the non-zero modes at 
one-loop making use of the formulae 

^ 1 _ 1 ^ 1 _ 1 

^ (2^ ~ 12 ' ^ ~ 4 ^ ^ ^ 

to sum over the Matsubara frequencies in the loop. This gives rise to the operators written 
in terms of zero modes as 

(Vi)2 )) = ^iV/?3/2tr {An? , ^ ^(( (V^f )) = -^Np'/hr {Aof , (4.9) 

Os^^ii {Vsf )) = -^N/3'/hr {Anf , 0,^^{{ {V,f )) = ^Np'/hv {A,f , (4.10) 

O5 = (( V5 )) = -^Np^/hr (Aof , Oe^ii Ve )) = -^Np^/hr (Aif . (4.11) 



^If we left the 't Hooft coupling constant A arbitrary, we would find that the expansion parameter is 
given by e^/A = a/A/T3. This is what one might have deduced on dimensional grounds, since the 't Hooft 
coupling constant has the dimension of (mass)^ in the present models. 
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Summing up these operators, we obtain 

= t.^j = -{^-l)NP'/'{iriA,f-iriAnf} . (4.12) 
Using this operator, we can evaluate the first term in eq. (4.7) as 

= X, r'/^ - - I) (X3 - X4) r-' + O(r-V') , (4.14) 

where the subscript "c" imphes that the connected part is taken, and we define the coeffi- 
cients X3 and Xi by 

X3 = (tr(i,)^.tr(i,-)^\ , X4 = (tr(A)2-tr(i^)2) . (4.15) 

\ / DR,c \ / DR,c 

The second term of eq. (4.7) can be calculated at the next-leading order using the propa- 
gator (4.3), and we get 

1 ^ (tr (x;xij) ^ 1 ^ ((tr (x;xij)) = ^T-' + o(r-5/2) _ ^^^g) 

Adding the two terms, we get the result at the next-leading order as 

{R^) = xiTV2 + - - I) (X3 - X4)} T-^ + 0(r-^/^) . (4.17) 

Let us calculate the internal energy (2.7) using (2.8). The operators and Si can be 
decomposed as 

^b = ^tr(F,.)2-^V6 + --- , ^f = -^V4 + ---, (4.18) 

where wc have omitted terms irrelevant at the next-leading order. The expectation values 
can be calculated as 

= ^X2 r + {l(d - i)xi - ^ - 1) (X5 - X6)} + o(r-) , (4.20) 
^'^^ " -2i^(<^(^^)'»)DR = -|^^^"^' + °(^") ' 

where we define the coefficients 

X5 = (tr(4)2.tr(ife)2\ , X6 = (tr(4)2-tr(i,5)') . (4.22) 

\ / DR,c \ / DR,c 
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Adding these terms, we get the result for the internal energy as 
Similarly the Polyakov line can be calculated as 

= 1 _ i. „T-3/^ + {^,, + ^(^^^ - I) - „)} T- + 0(T-"^) ,(4.24) 
where we define the coefficients X7 X8 by 



X7^(tr(iD)2-tr(io)2) 




(4.25) 



Thus we have obtained various quantities up to the next-leading order with the co- 
efficients Xij which can be obtained by Monte Carlo evaluation of the connected Green's 
functions in the DR model (3.2). In practice, we rewrite the Green's functions as described 
in Appendix B in order to increase the statistics. The values of Xi obtained in this way^ 
for d = 3,5,9 and for various arc summarized in table 1. In order to sec the large- AT 
behavior [25], we plot the values of Xi against 1/A^ for 12 < A < 32 in fig. 1. Wc observe 
that the data points for A = 16, 20, 32 lie on a straight line. This enables us to obtain the 
large- A extrapolated values shown in table 1. 

Using the coefficients Xi extrapolated to A" = oo, we evaluate the expressions (4.17), 
(4.23) and (4.24). The results for the bosonic case can be readily obtained by setting p = 
in the same expressions. In fig. 2 we show various quantities as a function of T for d = 3 
(left column) and d = 9 (right column) , respectively. The curves represent the results of the 
high temperature expansion^. The solid lines represent the leading order results, which are 
the same for the bosonic and supersymmetric cases. The dashed lines and the dash-dotted 
lines represent the next-leading order results for the bosonic case and the supersymmetric 
case, respectively. For comparison, we also plot the recent Monte Carlo data obtained at 
finite T for the bosonic model with d = 3 [23] and d = 9 [17], and for the supersymmetric 
model with d = 9 [15]. In both bosonic and supersymmetric cases, the high temperature 
expansion including the next-leading order terms seems to be valid at T > 2. 

^The heat-bath algorithm as described in ref. [25] has been used for Monte Carlo simulation. We have 
made 25M sweeps for rf = 3, 5 and 8M sweeps for d = 9 to obtain the data. We have checked that the 
expectation value appearing in the definition (3.6) of xi agrees with the previous results [25]. 

■'in refs. [15,17,23], we have presented the results of the high temperature expansion using the coefficients 
Xi obtained at the same A'^ as those used for Monte Carlo simulations at finite temperature. The quality of 
the agreement with the Monte Carlo data at high T is almost the same as in the present plots. 
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5. Summciry and discussions 

In this paper we have formulated the high temperature expansion for the supersymmet- 
ric matrix quantum mechanics with 4, 8 and 16 supercharges. While the non-zero modes 
become weakly coupled at high temperature, the zero modes remain strongly coupled and 
hence they have to be treated non-perturbatively. This makes the problem nontrivial, but 
we are able to obtain the next-leading order terms by evaluating connected Green's function 
in the bosonic IKKT model using Monte Carlo simulation. Since the fermions decouple at 
the leading order, it is highly motivated to carry out the next-leading order calculation. 
Indeed, our results including the next-leading order terms are in good agreement with the 
finite temperature calculations down to T ~ 2 in units of the 't Hooft coupling constant. 
Note also that Monte Carlo evaluation of the connected Green's functions in the bosonic 
IKKT model is by far easier than simulating the supersymmetric matrix quantum mechan- 
ics at finite temperature directly. This enables us to study the behavior at larger N and 
to make a reliable large-AT extrapolation. Our results confirm that the values of N used 
in finite temperature simulations are already large enough to probe the 't Hooft large-iV 
limit at high temperature. 

It is straightforward to extend our calculation to higher orders. For that, one needs 
to evaluate connected Green's functions with more than two operators inserted. That 
will require more statistics in Monte Carlo evaluation. For finite N, it is anticipated that 
the connected Green's functions with many insertions of the tr (A^)'^ operator would be 
divergent [27], and the order at which such divergence shows up would grow linearly with N. 
This property of the high temperature expansion is reminiscent of the infrared instability 
observed in Monte Carlo simulation of the supersymmetric model at finite temperature [15]. 
We note, however, that the divergence in the high temperature expansion occurs also in 
the bosonic case, in which the finite temperature calculations exhibit no such instability. 

It is worth while to generalize our formulation to higher dimensions. For instance, 
an interesting phase structure is expected in 2d U(A^) M = 8 super Yang-Mills theory 
on a finite torus. ^ In the strong coupUng and low temperature regime, the gauge/gravity 
duality predicts [10, 31] that there exists a phase transition corresponding to the black- 
string/black-hole transition [32] in the dual gravity theory. In the weak coupling and high 
temperature regime, on the other hand, one can study the theory by dimensionally reduced 
Id bosonic model. In rcf. [10], the phase transition observed in the latter regime has been 
conjectured to be a continuation of the phase transition in the former regime. In ref. [17] 
the critical region of the dimensionally reduced Id bosonic model has been studied more 
carefully, and a new phase characterized by the non-uniform eigenvalue distribution of 
the holonomy matrix has been discovered. Of particular interest from the viewpoint of 
the gauge/gravity duality is to investigate the fate of this new phase as one lowers the 

®The situation becomes more involved in dimensions higher than two. Recently cascade transitions 

from the black p-brane solution to the black (p — l)-brane solution have been found in the dual gravity 
theories [28]. On the other hand, the high temperature limit of the 4d U(iV) A/" = 4 super Yang-Mills theory 
on a finite torus, for instance, is described by the dimensionally reduced 3d bosonic model. Analogous 
cascade transitions were observed earlier in the large- iV pure Yang-Mills theory on a 3d torus [29] and in a 
related model [30]. 
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temperature. Calculations including the next-leading order terms in the high temperature 
expansion would be useful for such purposes. 
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Appendix 

A. Derivation of the formula for the internal energy 

In this Appendix we derive the formula (2.8), which is used to calculate the internal energy 
by the high temperature expansion. The case without fermions is given in ref. [17]. Let us 
first rewrite (2.7) as 

Z{/3)Af3-.o A/3 ' ^ ^ 

where (5' = f3 + A/3, and represent Z{f3') for later convenience as 

Z{I3') = J [VA'y [VX'y [V^%, e-^' , (A.2) 

where S' is obtained from S given in (2.1) by replacing f], t, A{t), Xi{t), ipa{t) with /?', t', 
A'{t'), X[{t'), '^a{t')- In order to relate Z{(3') to Z{l3), we consider the transformation 

t' = ^t, A'{t') = ^A{t), Xlit') = ^ X,{t) , ^p'{t')=m, (A.3) 

where the constant factors are motivated on dimensional grounds, and we have ['DX']^' = 
[VX]f3, [Dip'jp/ = \Pi^](} and [VA'jp/ = [DA] p. Under this transformation, the kinetic term 
in S' reduces to that in S, but the interaction term transforms non-trivially as 



'0 

This gives us the relation 



^ di'tr {[Xi{t'),X'^{t')]y = ^ dttr [[X,{t),XMy , (A.4) 

dttT{^^{t')[Xi'{t'),i;'f,{t')]) = (^^^ J^^dttr(^MMt),Mt)]) ■ (A.5) 



ZiP') = Zif3){l-N^Af3{S^ + £i) + 0{{Apf)) , (A.6) 

where the coefficients iSb and £f are defined by (2.9) and (2.10). Plugging these into (A.l), 
we get (2.8). Thus we are able to express the internal energy E in terms of the expectation 
values, which can be calculated directly by Monte Carlo simulation. 
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B. Increasing statistics by exploiting SO{D) symmetry 



In eqs. (3.6) and (3.7), we have rewritten the expectation values that define xi X2 
by exploiting the SO{D) symmetry of the DR model (3.2). Similar rewriting can be done 
also for the other coefficients Xi = 3, • • • ,8) defined in eqs. (4.15), (4.22) and (4.25) as 
presented below. In actual measurements in the Monte Carlo simulation, wc can increase 
the statistics considerably by using these forms instead of the original ones. 



fj.<iy 



«-p^EE(*'(^.^)^-''(^*)%,,e. 
« = ^E(''(^.)^-''(4)^)„,,e' 
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Table 1: The values of Xi (« = 1, 3, • • • , 8) for various d and N obtained by Monte Carlo simulation 
of the corresponding DR model. The values at A?" = oo are obtained by extrapolating the results 
for N = 16, 20, 32 as shown in fig. 1. 
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Figure 1: The coefficients Xi (i = 1, 3, • • • , 8) for d = 3,5,9 and 12 < TV < 32 evaluated by 
Monte Carlo simulation of the corresponding DR model are plotted against l/N"^. The straight 
lines represent fits to the expected large- A/' behavior a + b/N"^ using the A'' = 16, 20, 32 data points. 
The extrapolated values are shown in table 1 as results at A'' = oo. 
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Figure 2: Various quantities are shown for d = 3 (left column) and d = 9 (right column). The 
curves represent the results obtained by the high temperature expansion using the large- A'' extrap- 
olated values for Xi shown in table 1. The solid lines represent the leading order results, which are 
the same for the bosonic and supersymmetric cases. The dashed lines and the dash-dotted lines 
represent the next-leading order results for the bosonic case and the supersymmetric case, respec- 
tively. The circles and squares represent the results obtained by Monte Carlo simulation at finite 
T for the bosonic model with A'' = 16 [17,23] and for the supersymmetric model with A'' = 12 [15], 
respectively. 
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